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We discuss fits of cosmological dark energy models to the available data on high- 
redshift supernovae. We consider a conventional model with Cold Dark Matter and 
a cosmological constant (ACDM), a model invoking super-horizon perturbations 
(SHCDM) and models based on Liouville strings in which dark energy is provided 
by a rolling dilaton field (Q-cosmology). We show that a complete treatment of 
Q-cosmology requires a careful discussion of non-equilibrium situations (off-shell 
effects). The two main high-redshift supernova data sets give compatible constraints 
on ACDM and the other models. We recover the well-known result that ACDM fits 
very well the combined supernova data sets, as does the super-horizon model. We 
discuss the model-dependent off-shell corrections to the Q-cosmology model that are 
relevant to the supernova data, and show that this model fits the data equally well. 
This analysis could be extended to other aspects of cosmological phenomenology, 
in particular to the CMB and Baryon Acoustic Oscillations, which have so far been 
treated using on-shell models. 
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1 Introduction 

There is a plethora of astrophysical evidence today, from supernovae mea- 
surements, the spectrum of fluctuations in the cosmic microwave background 
[1-4] , baryon oscillations and other cosmological data, indicating that the ex- 
pansion of the Universe is currently accelerating. The energy budget of the 
Universe seems to be dominated at the present epoch by a mysterious dark en- 
ergy component, but the precise nature of this energy is still unknown. Many 
theoretical models provide possible explanations for the dark energy, rang- 
ing from a cosmological constant [5] to super-horizon perturbations [6] and 
time- varying quintessence scenarios [3], in which the dark energy is due to a 
smoothly varying (scalar) field which dominates cosmology in the present era, 
such as a time-dependent dilaton field [7]. 

The current astrophysical data are capable of placing severe constraints on 
the nature of the dark energy, whose equation of state may be determined by 
means of an appropriate global fit. Most of the analyses so far are based on 
effective four-dimensional Robertson- Walker Universes, which satisfy on-shell 
dynamical equations of motion of the Einstein- Friedman form. Even in modern 
approaches to brane cosmology, which are described by equations that deviate 
during early eras of the Universe from the standard Friedman equation (which 
is linear in the energy density), the underlying dynamics is assumed to be of 
classical equilibrium (on-shell) nature, in the sense that it satisfies a set of 
equations of motion derived from the appropriate minimization of an effective 
space-time Lagrangian. 

However, cosmology may not be an entirely classical equilibrium situation 
[7]. The initial Big Bang or other catastrophic cosmic event, which led to the 
initial rapid expansion of the Universe, may have caused a significant departure 
from classical equilibrium dynamics in the early Universe, whose signatures 
may still be present at later epochs including the present era. Three of us ( JE, 
NEM, DVN) have indeed proposed one specific model for the cosmological 
dark energy which is of this type, being associated with a rolling dilaton field 
that is a remnant of this non-equilibrium phase, described by a generic non- 
critical string theory [8-10]. We call this scenario 'Q-cosmology'. 

The central purpose of this paper is to confront the present cosmological 
data on high-redshift supernovae with non-equilibrium cosmologies [7-10], 
in which the dark energy relaxes at different rates, and compare the results 
with the predictions of the conventional ACDM model [5] . We also comment 
on a model with super-horizon perturbations superposed on an underlying 
Einstein- Friedman- Robertson- Walker Universe 1 [6]. 

As we explain in more detail below, care must be taken in interpreting the 



1 Our approach to this model is strictly phenomenological: we do address any of 
its deeper theoretical issues. 



Q-cosmology scenario proposed in Refs. [8-10]. Since such a non-equilibrium, 
non-classical theory is not described by the equations of motion derived by 
extremizing an effective space-time Lagrangian, one must use a more general 
formalism to make predictions that can be confronted with the current data. 
The approach we favour is formulated in the context of string/brane theory 
[11, 12], the best candidate theory of quantum gravity to date. Our approach 
is based on non-critical (Liouville) strings [7, 13, 14], which offer a mathemat- 
ically consistent way of incorporating time-dependent backgrounds in string 
theory. 

The basic idea behind such non-critical Liouville strings is the following. Usu- 
ally, in string perturbation theory, the target space dynamics is obtained from 
a stringy a- model [11] that describes the propagation of strings in classical 
target-space background fields, including the space-time metric itself. Con- 
sistency of the theory requires conformal invariance on the world sheet, in 
which case the target-space physics is independent of the scale characterising 
the underlying two-dimensional dynamics. These conformal invariance condi- 
tions lead to a set of target-space equations for the various background fields, 
which correspond to the Einstein/matter equations derived from an appropri- 
ate low-energy effective action that is invariant under general coordinate trans- 
formations. Unfortunately, one cannot incorporate in this way time-dependent 
cosmological backgrounds in string theory, since, to low orders in a perturba- 
tive expansion in the Regge slope a', the conformal invariance condition for 
the metric field would require a Ricci-flat target-space manifold, whereas a 
cosmological background necessarily has a non- vanishing Ricci tensor. 

To remedy this defect, and thus be able to describe a time-dependent cos- 
mological background in string theory, the authors of Ref. [7] suggested that 
a non-trivial role should be played by a time-dependent dilaton background. 
This approach leads to strings living in numbers of dimensions different from 
the customary critical number, and was in fact the first physical application of 
non-critical strings [13]. The approach of Ref. [7] was subsequently extended 
[8-10, 14] to incorporate off-shell quantum effects and non-conformal string 
backgrounds describing other non-equilibrium cosmological situations, includ- 
ing catastrophic cosmic events such as the collision of two brane worlds, etc. 

In our discussion of such models in this paper, we first review briefly their 
main predictions. We then demonstrate the importance of off-shell corrections 
to the naive equations of motion of non-critical Liouville string cosmology, 
which cause us to modify the naive equations of motion of Q-cosmology. We 
then demonstrate that the available supernova data are compatible with such 
non-critical-string-based cosmologies [8, 10]. 

The structure of the article is as follows. In Section 2 we summarise the basic 
equations of the first set of models we study, namely ACDM and the super- 
horizon model. Then, in Section 3 we discuss the Liouville Q-cosmology model 



of Refs. [8-10] . In particular, we discuss technical aspects of Q-cosmology, plac- 
ing the emphasis on the off-shell corrections and deriving the fitting formulae 
to be used in our analysis. This part may be omitted by the reader who is 
more interested in the confrontation of this and other models with the avail- 
able high-redshift supernova data that we present in Section 4. We find that 
the ACDM model fits the data very well, as does the super-horizon model. We 
also show that the off-shell corrections to the Q-cosmology enable this model 
to be reconciled with the supernova data on equal footing with ACDM. 2 Our 
results and the possible outlook are summarised in Section 5. 

2 Basic Formulae for the Phenomenological Analysis of On-Shell 
Dark Energy Models 

In this Section we review the basic formulae needed in the phenomenological 
analysis of observational constraints on the cosmological parameters derived 
from supernova data and (in the case of the ACDM) also baryon acoustic 
oscillations, in the context of on-shell models (i.e., excluding the Q-cosmology 
model to be discussed later). The relevant dynamics of any such cosmological 
Robertson- Walker model may be lumped into the expression for the Hubble 
parameter H(z) in terms of the various components of the energy density 
of the Universe, represented in units of the (present-day) critical density for 
a spatially-flat Universe, Qi = p®/p®, where % — M for matter, including 
dark matter, % = A for constant dark energy, i.e., a cosmological constant 
A, etc. Later in this work we discuss fits to the on-shell (ACDM and super- 
horizon) and off-shell (Q-cosmology) models discussed in the current and the 
next Sections, respectively. 

2. 1 Standard A CDM Model [5] 

In this case, we consider a spatially flat Robertson- Walker Universe containing 
dust-like matter that can be described as an ideal fluid, and a cosmological 
constant A > corresponding to a constant equation of state of the form: 

P = wop , (1) 

where p(p) is the pressure (energy density) and wq = — 1. The standard 
Einstein-Friedmann equations for such a Universe yield: 

H(z) = H fa M (l + zf + n A (l + zf^A 1 ' 2 , (2) 



where Hq is the present-day value of the Hubble constant. 



2 A more complete comparison of this model with data on the cosmic microwave 
background (CMB) and baryon acoustic oscillations (BAO) would require further 
formal developments that lie beyond the scope of this paper. 



2.2 Super-Horizon Perturbation Model [6] 

In this model, the Universe is assumed to be filled with only non-relativistic 
matter, and there is no Dark Energy of any sort. The metric takes the form: 

ds 2 = -dt 2 + a 2 (t)e- 2 ^^ t) S lJ dx i dx j , (3) 

where \l/(x, t) denotes the gravitational potential. The authors of Ref. [6] 
considered perturbations ty(x,t) that correspond to conformal contraction or 
stretching of the three-dimensional Euclidean space, {x\ i — 1,2, 3}. Splitting 
the potential \I/ = tyg + \I/ S into contributions from long-wavelength (super- 
horizon) modes, \J/^, and short- wavelength (sub-horizon) modes, \l/ s , and ig- 
noring the latter within our Hubble volume, one obtains [6]: 

ds 2 ~ -dt 2 + a 2 {t)5 lJ dx i dx j , a(t) = a(t)e-* dt)+ * m , (4) 

with the normalization a — Oo — 1- 

The analysis of Ref. [6] used the appropriate Einstein's equations to demon- 
strate that, up to trivial rescaling of spatial coordinates and absorbing the 
time-independent part of \t>£, the time dependence of ^g(t) may be cast in the 
form: 

**(*) ~ a(t)V m . (5) 

One therefore finds that the expansion rate of the Universe, as observed by a 
local observer who is restricted to observe within our Hubble volume, is given 
by: _ 

ff w-- l!? f-r^(«- 3/2 -«- 1/2 *»)- < 6 » 

where (1 + z)~ x = a(t). Combining (4) and (5), one may extract the following 
relation between a and z: 

l + z = a' 1 e^ a ' 1 ^ m . (7) 

This model predicts acceleration of the Universe at the current era: 

3/2 - aVgo/2 



q = -l + 
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It can easily be shown by direct comparison with conventional models that 
the super-horizon perturbations mimic a dark energy model with equation of 

state 

Wde^de = g U ~ ^ J ■ ( 9 ) 

In order to evaluate the constraints imposed by astrophysical data on this 
super-horizon model, we use the following integral of the Hubble parameter 



E(z) = H(z)/H over z: 

JO £/(;z) ia(z) 

v^ e-^°^E§|(Erf(y=^o~) - Erf(^-tf«,a(z))), if *«, < 
2(l- A /^)Y if^ = (10) 



v^ e"**^ Erfi(v^) - Erfi(*Moa(z)) , if ^o > 



where Erf and Erfi are the standard and imaginary error function, respec- 
tively. We assume values |^o| < 1, in order for perturbation theory to remain 
valid [6]. We reiterate that we do not address any of the deeper theoretical 
issues of this model. 

2.3 A Note on Relaxing Dark Energy Models 

There are many dark energy models in which the energy density relaxes in 
time. For example, relaxation at a rate proportional to 1/t 2 , where t is the 
cosmic time, may arise in certain non-critical string theory models [14] . In such 
a model, the role of the quintessence field is played by the dilaton field [7]. 
However, at least in the string theory context, a complete treatment of such a 
rolling dilaton field requires a non-trivial discussion of off-shell physics, which 
we present in the next section. For the moment, we neglect this complication 
(as might be valid for some non-string models of relaxing dark energy), and 
assume a naive on-shell physical framework in which the dark energy density 
relaxes according to a generic power law l/t n . Within this general framework, 
n = corresponds to a cosmological constant A and n = 2 corresponds to a 
naive interpretation of the Q-cosmology model of Ref. [8], in which only the 
conventional type of matter exists and the off-shell terms are not taken into 
account. Such a relaxing dark energy model would predict a Hubble parameter 
of the form: 

H(z) = H (fi M (l + zf + %{l + z) n ) 1 ' 2 , (11) 

where Q,p represents the present-day dilaton dark energy in the on-shell for- 
mulation. 3 In such a model with n = 2, the rolling dilaton dark energy would 
resemble a spatial curvature contribution, and such naive on-shell models are 
easily excluded on the same observational grounds that support experimen- 
tally the spatial flatness of the Universe. However, n = is certainly allowed, 
and it is important to know the permissible range of n. There is a formal corre- 
spondence between (11) and a generalized ACDM model in which the vacuum 
energy A obeys an equation of state with w ^ —1 (2), namely n = 3(1 + w ). 
However, the underlying physics may be entirely different. For models with 



3 However, we emphasize again that this on-shell treatment is incorrect in the 
context of non-critical string theory, as we discuss below. 



relaxing dark energy (11), the standard data analysis shows that both the 
'gold' and the SNLS dataset favour a value of the parameter n appearing that 
is <C 2. Indeed, the combined data prefer n ~ (w ~ —1), i.e., a (nearly) 
constant vacuum energy, as in the ACDM model. 

3 Non-Critical Liouville String Q-Cosmologies 

In this Section we explore in more detail Q-cosmologies with relaxing vacuum 
energy in the context of Liouville strings, which provide a consistent treatment 
of off-shell, out-of-equilibrium effects. For the benefit of the reader, we review 
here the basic formalism and summarize the results, giving details of their 
derivation in the next two subsections, that may be skipped by a reader not 
interested in technical aspects. There are many specific cosmological models 
that may be formulated within the general framework of Liouville strings: 
see [8-10] for an extended discussion and references. Here we limit ourselves to 
the main features of this cosmological framework, which do not depend on the 
details of the underlying microscopic string/brane models or compactification. 
We simply assume the theory to have been compactified (somehow) to a four- 
dimensional Robert son- Walker space-time, with a scale factor a(t^) in the 
so-called Einstein cosmic frame [7]. 

3.1 General Formalism of Off-Shell Liouville Q- Cosmologies 

In the absence of matter, the Liouville-dressing approach of Ref. [14] and the 
(dynamical) identification of the Liouville zero mode with target time lead to 
generalized conformal invariance conditions for the fields of the gravitational 
multiplet of the string propagating in a four-dimensional background: 

g' + Qg^-F, (12) 

where the overdot denotes differentiation with respect to the world-sheet zero 
mode of the Liouville field 0, the set of g l = {G^, 0, . . . } are string back- 
ground fields (metric, dilaton, etc.), the f3 l are the respective world-sheet 
renormalization-group /5-functions, and Q is the square root of the central- 
charge deficit that describes the departure from criticality. The negative sign 
in front of f3 l on the right-hand side of (12) corresponds to supercriticality for 
the string, in which case the central charge exceed its critical (conformally- 
invariant) value. It is in this case that the Liouville mode has a time-like 
signature and can be identified with the target time [7, 14]. 

In critical strings, the vanishing of the f3 l would correspond to the ordinary 
Einstein/matter equations of the low-energy field theory derived from strings. 
The vanishing conditions correspond to an equilibrium situation. In non- 
critical string, however, these equations are replaced by (12), which express the 
restoration of conformal invariance by the Liouville mode. The physical rea- 
sons for an initial departure from conformal invariance might be a catastrophic 



cosmic event, such as the collision of two brane worlds or a Big Bang [8]. Such 
an event would cause a departure from equilibrium, which may not be inappro- 
priate for the early era of the Universe. Following the identification of target 
time with the Liouville mode, the generalized conformal invariance equations 
(12) yield new dynamical equations that describe the non-equilibrium dynam- 
ics of the Liouville strings [8-10, 14]. 

If one further assumes that, after compactification to four space-time dimen- 
sions, one obtains a spatially-flat cosmological Robertson- Walker metric back- 
ground, e.g., in the context of a brane model, the equations (12) define the 
dynamics of the resulting 'Liouville Q-cosmology', as we term such a non- 
critical string cosmology. The equations (12) are not the ordinary equations 
of motion corresponding to a four- dimensional gravitational effective action, 
but describe the dynamics of an off-shell relaxation process. 

As explained in Ref. [8], the inclusion of matter does not change qualitatively 
this non-equilibrium, off-shell nature of such a model. Apart from the formal 
addition of extra terms in the gravitational /9-functions, resulting from the 
coupling of matter to the gravitational field, which are needed for general 
coordinate invariance in the target-space, in most models of physical interest 
may be included in such a way that the main source of non-criticality is the 
gravitational and moduli sector, while the matter f3 l functions themselves are 
(almost) vanishing. In such a case, the matter fields feel the presence of the 
off-shell Liouville terms on the left-hand-side of (12) via their gravitational 
couplings. 

The presence of these off-shell Liouville terms, which are exclusive to non- 
critical strings, results in a modification of the critical energy density condition 
for the Liouville cosmology, which ensures spatial flatness. This critical density 
is therefore different from that in conventional (on-shell) Friedman-Robertson- 
Walker (FRW) cosmologies. 

The essential formalism is that of Ref. [7], in which all physically relevant 
quantities should be expressed in the Einstein frame, in terms of the Einstein 
cosmic time. The four-dimensional matter action (including radiation fields) 
couples to the dilaton field non-trivially, in a way that is specific to the various 
matter species, as a result of purely stringy properties of the effective action 
[11]. A generic a-model-frame effective four-dimensional action with dilaton 
potential V ((/>), which could even include higher string-loop corrections, has 
the form: 

5 (4) = — [ d 4 xV^G[e-* w R(G) + Z(0)(V0) 2 + 2a'V((j)) . . . 1- 

1 f , 1 ( 13 ) 

— J d 4 xVG——F^ - I m ((f), G, matter) , 

in the notation of Ref. [15], with the various factors \l/, Z, a encoding informa- 



tion about higher string loop corrections. Also, F^ v denotes the radiation field 
strength and I m ((f>, G, matter) represents matter contributions, which couple to 
the dilaton in a manner dictated by string theory, with specific scaling laws 
[11] under shifts of the dilaton field — » + const. In cases [8] where only the 
string tree level plays a role for late times, the various form factors simplify, 
e.g., ^(0) = 20, Z(<f)) = 4, etc. However, for reasons of generality, here we 
keep the form (13). As we shall see, there is sensitivity in the present data to 
terms in the dilaton potential with higher powers of the string coupling. When 
higher loop corrections are important, these factors have a complicated form, 
for instance one has e*^ = CQe~ 2 ^ + Ci + C2e 2 ^ + . . . , with various constants q, 
with the powers of the square of the string coupling g 2 = e 2 ^ counting closed 
string loops, as appropriate for the gravitational multiplet. For simplicity, in 
the present work we ignore the antisymmetric-tensor four- dimensional field, 
which, as discussed in Ref. [7], corresponds to an axion field. 

Higher loop corrections may result in a modified dilaton potential for the 
noncritical string theory of the form 

V = 2Q 2 e 2 * + V , (14) 

where the central charge deficit Q 2 encodes the microscopic non-equilibrium 
physics [8, 9]. In the models we are considering in this work, Q 2 > and the 
string theory is supercritical [7]. 

Assuming a normal fluid form for matter or radiation, with stress tensor T vE = 
diag \—p,p$i) in the Einstein frame [7], we obtain the following gravitational 
equations of motion in our case (we work here in units Mp = 1/8ttGn = 1, 
where Mp is the four-dimensional Planck constant) [10]: 







dn 2/ 
2-t— = -p m ~ P<p - Pm - P<t> + a [t E )Jii , i ■■ 
diE 


= 1,2, 


»=H 2( £ ),+ H' * s H 2( £ )2 - 


-y(0) 


rf 2 , ou d<P i 1 dV , 1 f N 3J tl 
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(15) 



J4, 



where t E is the Einstein-frame Robertson- Walker cosmic time [7] and the 
potential of the dilaton assumes the form (14), including loop corrections V 
(14). We note that the equations (15) differ from those holding in on-shell 
cosmologies by the Liouville out-of-equilibrium contributions J7, which are 
exclusive to our treatment [8, 10, 14]. 

The equations (15) need to be supplemented by the appropriate Curci-Paffuti 
equation [16] expressing the renormalisability of the cr-model fields [9, 10], 



which is also exclusive to our non-critical string approach. This equation re- 
lates the dilaton world-sheet /5-function to the rest of the /9-functions in the 
problem (for the graviton, etc.), and essentially expresses the (time-dependent) 
central charge deficit in terms of the physical fields in the problem, 



^ = -6e- 2 *(if + 0)4 
aiE a 



-6e- 2 *(# + 0)^ + ... , (16) 



where the overdot henceforward denotes a derivative with respect to the cosmic 
time i_e, the '. . . ' indicate terms due to higher order in string loops, and the 
values of J$ and Ju are given by 



J = e- 2 ^(0-0 2 + Qe < V), 

Ju = 2a 2 U + 3#0 + 2 + (1 - q)H 2 + Qe*(0 + H)) , i = 1, 2, 3, 



(17) 



with q the deceleration parameter: 



. , aa dH/dtE (1 + z)dH/dz 

*--> B -? = -7?4f - ' - ww - 1 - (18) 

For completeness, we mention that the equations (15) lead, after standard 
manipulations [8, 10], to the (non-)conservation equation of matter in the 
presence of the non-equilibrium contributions: 

^m , Qff , . x ldQdV((f>) d(f) R/W, #\-2/r f. a \ 

-^ + mp m + Pm ) + - 2 ^- E -^q—^ (p m -3 Pm )=6(H+— )a J, (19) 

These equations are used below to obtain a general expression for the Hubble 
parameter as a function of the redshift, which we use later in our fit to data. 



3.2 Derivation of Basic Phenomenological Formulae for Q- Cosmologies 

The above analysis contains non-linear equations that are hard to solve an- 
alytically. A detailed numerical analysis is performed in Ref. [10], where we 
refer the interested reader for details. Instead, for our phenomenological pur- 
poses here, we shall attempt an approximate analytic solution at late epochs 
of the Universe, appropriate for our phenomenological fits. Throughout this 
work we work in units where the present scale factor is normalised to one, 
a = 1. 

To solve (19) in the various epochs of the Universe, it is convenient first to split 
the energy density of matter into radiation p r , baryonic pb and dark-matter 
Pd components, and to use simple equations of state for the dilaton fluid, as 



10 



supported by our theoretical model [8-10]: 

Pm = Pr + Pb + Pd = Pr + PM , 
Pb = 0, p d = W d p d , p r = -p r , p^ = W^p^ . (20) 

We have assumed a non-trivial equation of state Wd for the dark matter com- 
ponent, which as we discuss below, stems from the fact that its scaling with 
the scale factor a is exotic, different from the ordinary dust scaling a~ 3 . One 
can split the matter evolution equation (19) into various components, using 
(20). 

In the solution of Ref. [9] without matter, at late eras the dilaton <j) and the 
scale factor a(t) vary as follows with the cosmic time £#: 

(j) = -lna(t E ) , 
o = a 1 (l+ 7 2 4) 1/2 , (21) 

where a±, 7 are appropriate positive constants. 4 The form of the scale factor 
(21) implies the following form for the deceleration parameter at late epochs 
of the Universe: 

q(z) = -l/( 7 2 t|) , (22) 

The values of the constants are such that the present era in the history of 
the Universe is characterised by 7^ = 0(1), as follows from the fact that 
the deceleration parameter of the Universe is currently observed [2b] to be 
go = -0.61. 

We notice that the asymptotic behaviour of the model of Ref. [9] without 
matter, at very late times, indicates that 

+ #~O, (23) 

and thus 

H ~ a -1 , (24) 

implying that q ~ 1/aH, varying very little at late epochs. 

In the present case, where matter effects are important, the formulae need 
some modification. In general, the dilaton and the scale factor differ from the 
particular form (21) [10], in which case one obtains a much more complicated 
behaviour for q(z) than the one given in (22), which can only be studied 
numerically. A complete numerical analysis is performed in Ref. [10], where 
we refer the reader for details. We recall, for completeness, that the analysis of 



4 In the model of Ref. [9] they are related to the flux in the extra compactifled 
dimensions of the appropriate string theory, but they are viewed here as arbitrary 
constants in our more general setting. 
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Ref. [10] indicates evidence for a past deceleration of the Universe at redshifts 
larger than z* = 0.37. 

For our generic phenomenological purposes in this work we can make several 
physically sensible simplifying assumptions, which allow us to carry some of 
the results of Ref. [9] through to the present case, allowing for some analytic 
treatment. Specifically, the dilaton and scale factors are assumed to take the 
form 

= -\na(t E ) + 0(g 2 s ) , (25) 

as dictated by a general analysis of the dilaton equation in the model. 5 Notice 
that, in view of the general discussion in Ref. [10], we no longer maintain the 

1 /9 

explicit form of a — a\ (1 + l 2 t%) used in the matter- less model of Ref. [9], 
however we do maintain the scaling relation (24), which, according to the 
numerical analysis of Ref. [10] , seems to be a good approximation in the eras 
of interest to us here. 

In general [10], there may be a more complicated behaviour than (23) or (24). 
For instance, loop and other corrections may yield power-series corrections in 
the (perturbative) string coupling 

<P + H = 0(Hg 2 s ) + ... , g a = e+ . (26) 

Such corrections could lead to appreciable deviations in the behaviours of var- 
ious fields and other cosmological parameters even at the present era. In fact, 
the numerical analysis of Ref. [10] indicates evidence for a past deceleration of 
the Universe at redshifts smaller than z* = 0.37. Nevertheless, even in these 
more general models, the Hubble parameter maintains [10] a scaling similar 
to (24) at the late eras of interest to us here, which implies that some of the 
qualitative features of our analytic approach used here are valid, to a good 
approximation, in these more complicated cases. 

On using (15) and performing some elementary manipulations, we easily ob- 
tain: 

p M ~ 2(2g - 1)H 2 ~\^ + \v + \e^J, . (27) 

A consistent assumption, as we shall demonstrate now, is that the matter 
energy density, pm, includes dust (baryonic matter) and exotic scaling (dark) 
components, which can behave like radiation-like contributions: 

Pm ~ Pdust a ~ 3 + Pcxotic a_4 + • • • , (28) 

Indeed, taking into account (24), we observe that the above assumption is 
consistent with the conservation equation (19) provided that the central charge 



We thank V. Georgalas for an informative discussion on this point. 
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deficit at late eras scales as 

Q 2 {a) ~ Q2 + Pdust (2g) 

a 

with Ql > a constant, namely the central charge of the asymptotic conformal 
field theory to which the string model flows at infinite time [8, 9]. 

Notice that, if one allowed an exotic a~ 2 scaling of matter of the form p^_ like a~ 2 , 
which characterises the dilaton dark energy p^, as we discuss below, then Q 2 
should contain also terms of the form p9_ like ln(a 2 ). This would lead, in gen- 
eral, to logarithmic scaling a _2 ln(a) of the various energy density components, 
which would complicate the situation, as far as the validity of the approxima- 
tion (24) is concerned. Although such terms may exist in some models, we do 
not consider them in this work. Here we assume the scaling (28) for the matter 
energy densities in the eras of interest to us. For more complete, numerical 
treatments we refer the interested reader to Ref. [10]. 

To proceed with a consistency check and estimates on the scaling behaviour 
of p m , P(f>, it is necessary to make some assumptions about the form of the 
potential V. As already mentioned, we assume that it is generated by higher 
string loops. In brane-inspired models such as those considered in Refs. [8, 9], 
one has open strings living on the brane world. The corrections then will be 
proportional to an extra power of the string coupling g s = e*. Assuming loop 
corrections of the form 6 

V -ae^ + fie^ + ... , a, (3, ... = const., (30) 

and taking into account (29), we observe that the dark matter density (27) 
yields: 



2(2g -1)+Q 2 pLt P 1 | 1« 
a 2 a 3 2 a 4 2 



Pm ~ 5 1 s — -7 + -e ^J^ + ... . (61) 



As we see from (26), (24), (17), (29), the quantity J^ = 0(a~ 2 ) is negative. We 
also note that the late-era scaling (28), which is required for consistency (on 
account of (29)), can easily be guaranteed, provided the terms 2(2g — 1) + Q 2 
which otherwise would yield a a~ 2 scaling are suppressed for the redshift region 
< z < 2 of interest to us. This would require a Ql of order 2|(2g — 1)1) 
according to our assumption that q does not change much in the region in 
which we are interested. 7 



6 The dots indicate higher-order corrections in the string coupling and also possible 
dilaton- independent contributions . 

7 We remark at this point that the actual data may deviate from this for redshifts 
less than 0.5, indicating past deceleration of the Universe. However, even in such a 
case, the absence of a~ 2 terms in matter today might be a consistent assumption. 
In fact, the numerical solution of Ref. [10] indicates that the matter does develop 
an exotic scaling cT 2 but at much later eras than the present one. 
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Notice that the dust-like scaling does not receive contributions from the string 
loop corrections, and in fact a comparison between (31) and (27) yields a 
mathematically consistent result with p^ust f ree parameter to be constrained 
by fitting data and/or other microscopic model considerations. 

In contrast, the exotic scaling a -4 receives contributions from the loop cor- 
rections, proportional to /3, whose sign cannot be fixed by our generic consid- 
erations, and depends on the details of the microscopic string/brane theory 
model, and contributions due to the off-shell Liouville terms e 2< ^j7^ which are 
positive according to our assumption that Q* > 0: 



o _P _ \q\ + Q* 

exotic J o o 



(32) 



The asymptotic behaviour of the dilaton dark energy may be evaluated from 
the definition of p^ (15): 



12 V tt2 Ql Pdust + a/2 P 

2 a 1 a 6 2 a 4 

0(a- 2 ) + 0(a- 3 ) + ^ a' 4 + . . . , 



(33) 



where the dots have the same meaning as in (30). It is interesting to observe 
that the dilaton dark energy component contains a part scaling like a~ 2 , which 
is positive and of order (1 + Ql)a~ 2 , but also parts scaling like a~ 3 and a~ 4 , 
which appear also in the matter energy density (27). However, in contrast to 
the matter case (31), the 'dust '-like contributions to the dilaton dark energy 
do not have a fixed sign or magnitude, as they depend on the string loop 
correction parameter a, 



J) 



b,3 



2p° us t + «/2~(l + «)/2. (34) 



When combined with the corresponding parts in Pm, in order to obtain the 
precise scaling of the Hubble parameter with the redshift z, H(z), from the 
first of equations (15), we observe that 

p, + p M ^ \0{a' 2 )\ + l4U±g a -3 + (^-\Sl±9l + . . . ) a~\ (35) 

where we have been careful to indicate the relative signs of the various terms, 
as predicted by the model with the approximations made so far. 8 Notice the 
independence of the exotic scaling on the loop parameter /3, which does not 
enter the expression for H(z). Thus, it is a feature of this Liouville model that 



8 One may have deviations from (26) [10], resulting in extra contributions to the 
a -4 scaling, denoted by dots. Thus, the exotic (dark) matter Q$ could be positive. 
This case seems to be favoured by the data, as discussed in the text, and it is in 
agreement with the numerical analysis of Ref. [10]. 
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the sign of the dust-like contribution appearing in H(z) is not fixed, since such 
terms receive contributions from the dilaton dark energy. Indeed, in modern 
string/membrane models, one may obtain negative dust-like contributions to 
the dilaton potential from compactification of Kaluza-Klein graviton modes in 
brane models, etc. [17]. In contrast, the exotic scaling terms appear positive. 

When fitting the data, in order to distinguish ordinary matter from dilaton 
dark energy effects that scale similarly with the redshift, one would have to 
perform also measurements of the equation of state of each component, which 
we do not attempt in this work. In the specific model [9] considered here as 
our pilot study, the values of the various fij appearing in the expression for 
H(z) are in principle predicted by the underlying microscopic dynamics of the 
model, as we have seen above, but we leave this exercise for future work. 

However, since some of the Qj involve highly model-dependent coefficients of 
string loop corrections and other, currently unknown, details of microscopic 
models, and, moreover, in the above derivation several assumptions have been 
made, e.g., about the form of the potential for the dilaton, the form of the 
dilaton and the scale factor at late epochs (21), etc., which may not always 
be valid [10], we consider at present the various parameters appearing in the 
expression for H(z) as arbitrary parameters to be fixed by the data. To be 
even more general, we keep the exponent of the exotic scaling as a free fitting 
parameter as well. From the above considerations, in particular (28), (33), it 
is also clear that in general f2d us t may be negative, as its sign depends on the 
sign of the corrections (30) to the dilaton potential. 

The final result of our parametrisation for H(z) in the Q-cosmology frame- 
work, which will be compared with the astrophysical data in the next section, 
in order to obtain information on the basic cosmological parameters of the 
model, is therefore: 

H(z) = H {n 3 (i + z) 3 + n s (i + z) s + n 2 (i + z) 2 ) 1/2 , (36) 

with the densities ^2,3,s corresponding to present-day values (z = 0) and 

n 3 + n s + n 2 = 1 , (37) 

The reader should recall that in the model described in this section 5 = 4. 
However, in view of the various approximations employed in our analytic treat- 
ment of Q-cosmology above, which, as already mentioned, may not be valid 
in the present era, when matter contributions are important, the exponent S 
is treated from now on as a fitting parameter. The above formulae are valid 
for late eras, such as the ones pertinent to the supernova and other data 
(0 < z < 2) that we use in this work. 

We stress once again that the various Qi contain contributions from both dark 
energy and matter energy densities. As explained previously, Q 3 does not 
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merely represent ordinary matter effects, but also receives contributions from 
the dilaton dark energy. In fact, the sign of f2 3 depends on details of the under- 
lying theory, and it could even be negative. For instance, Kaluza-Klein graviton 
modes in certain brane-inspired models [17] yield negative dust contributions. 
In a similar vein, the exotic contributions scaling as (1 + z) s are affected by 
the off-shell Liouville terms of Q-cosmology. It is because of the similar scaling 
behaviours of dark matter and dilaton dark energy that we reverted to the 
notation Qi, i = 2, 3, 5 in (36). To disentangle the ordinary matter and dilaton 
contributions one may have to resort to further studies on the equation of 
state of the various components, which we do not study in this article. More 
generally, one could have included a cosmological constant fi^ contribution in 
(36), which may be induced in certain brane- world inspired models. We do 
not do so in this work, as our primary interest is to fit Q-cosmology models 
[8-10], which are characterised by dark energy densities that relax to zero. 



4 Data Analysis 

4-1 Astrophysical Observables 

We now discuss the observables available to test the above models using mainly 
data on high-redshift supernovae and (to a limited extent) baryon acoustic 
oscillations. 

We use the supernovae data reported in Refs. [4, 18], which are given in terms 
of the distance modulus 



/x = 51ogd L + 25 . (3£ 



We use for d^ the standard formula relating the luminosity distance to the 
redshift z 

dz' 

/o Hiz 7 



d L = c(l + z) /_ — -, (39) 



where H is given by (2), (6) and (36) in the three models discussed in the pre- 
vious sections, namely ACDM, super-horizon and Q-cosmology respectively. 
We note that this observable depends on the expansion history of the Universe 
from z to the present epoch, and recall that most of the available supernovae 
have z < 1, though there is a handful with larger values of z. 

For standard on-shell cosmologies, complementary information is provided by 
the data on baryon acoustic oscillations (BAO) [19, 20], which show up in 
the galaxy-galaxy correlation function at z ~ 0.35. However, for reasons to 
be discussed later on in the article, the application of such an analysis to the 
off-shell Q-cosmology model is an open issue, since the underlying theoretical 
framework needs to be re-evaluated. 

The formulae for H(z) in the various models discussed above, namely (2), (6) 
and (36), serve as the basis for our analysis below. 
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4-2 Experimental Fits 

We now present our fits of the above models to the presently available astro- 
physical data on high-redshift supernovae [4, 18], and (in the case of ACDM) 
also baryon acoustic oscillations (BAO) [19, 20]. We start our discussion with 
the ACDM and Super-Hubble-Horizon models, and then proceed to the case 
of off-shell Liouville Q-cosmology models [8-10]. 

Measurements [1, 4] are available of the distance moduli (38) for 157 super- 
novae in a so-called 'gold' sample [4], observed by ground- based facilities and 
the Hubble Space Telescope: see Fig. 1. We compare these with the measure- 
ments of 71 other high-redshift supernovae published more recently [18] by 
the Supernova Legacy Survey (SNLS), which are accompanied by a reference 
sample of 44 nearby SN la, yielding a total of 115 data points: see Fig. 2. For 
both samples, the data are expressed in terms of the observed distance mod- 
ulus n b derived from the SN la light curves and the corresponding redshift z. 
The best-fit cosmological parameters are then extracted by minimizing: 

2 _ >T^ K.f'ob^ — t l th\ z i)) (AV\\ 

X - 2^ 2,2 ' y^ u ) 

i u ob,i * u int 

where /j, j, : i and a ^i are the measured value for the distance modulus and 
the corresponding error for a specific supernova i, respectively, Oi n t is the 
intrinsic dispersion of the absolute magnitudes of the supernovae, and Hth{zi) 
is the theoretical prediction for a given model as calculated by (38), (39), 
where d^ is given in units of megaparsecs. This fit and all subsequent analyses 
are performed with the ROOT [21] implementation of the Minuit function 
minimization and error analysis code [22]. 

Both the data and the predictions of the various models are expressed in the 
following as residuals, A/i, from their predicted values in an empty Universe 
(^M = 0): 

d ( r ty) = £-z(2 + z). (41) 

In addition to presenting the 'gold' and the SNLS data sets respectively, Fig- 
ures 1 and 2 also display the predictions of various cosmological models to be 
detailed later on in the article. 

In addition to the 157 SN la belonging to the 'gold' sample analysed in 
Ref. [4], a so-called 'silver' dataset of 29 SN la is also available, with lesser 
spectrometric and photometric quality compared to the 'gold' one. Although 
the results presented here were obtained with the 'gold' dataset, the analy- 
sis was also repeated including the 'silver' supernovae data — with a total of 
186 supernovae — with results comparable to those presented here, proving the 
robustness of the analysis. The analysis was also repeated for the combined 
sample of the 'gold' (157 supernovae) plus the SNLS sample (71 supernovae), 
which yields a sample of 228 supernovae. 
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Fig. 1. Residual magnitude of the 'gold' dataset supernovae as a function of the red- 
shift for the best-fit values of the cosmological parameters. The following models are 
shown: (i) Empty Universe (black solid line); (ii) Universe with matter only, Qm = 1 
(black dashed-dotted line); (iii) ACDM model for Qm = 0.287 (red long-dashed 
line); (iv) the super-horizon model for ^iq = —0.79 (blue short-dashed line); and 
(v) off-shell Q-cosmology model for $^3 = —3.6 and O4 = 1.2 (green thick solid line). 
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Fig. 2. Residual magnitude of the SNLS dataset supernovae as a function of the red- 
shift for the best-fit values of the cosmological parameters. The following models are 
shown: (i) Empty Universe (black solid line); (ii) Universe with matter only Qyi = 1 
(black dashed-dotted line); (iii) ACDM model for f^M = 0.265 (red long-dashed 
line); (iv) the super-horizon model for ^iq = —0.94 (blue short-dashed line); and 
(v) off-shell Q-cosmology model for ^3 = —3.0 and 0,4 = 1.0 (green thick solid line). 



The x 2 values obtained for the best fits to the various models are discussed in 
the remainder of this Section. 
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Fig. 3. Confidence-level contour plots for the ACDM model, (i) Left panel: 'gold' 
sample of 157 supernovae [4] and SNLS sample of 71+44 supernovae [18], (ii) right 
panel: combined 'gold' and SNLS sample of 157+71 supernovae. The 68.3% C.L. 
region extracted from the BAO analysis is also superimposed. 

4-2.1 Cosmological- constant model 

For the ACDM model [5], assuming a flat Universe, we find the \ 2 and pa- 
rameter values listed in Table 1, which are in agreement with those presented 
in Ref. [4, 18]. The corresponding confidence limits in the (Om, ^a) plane are 
shown in Fig. 3. 



Table 1 

Fits to the ACDM model parameter Qm, assuming a flat Universe. We compare the 
values favored by the 'gold' and the SNLS data sets, and give in the third row the 
results of a fit to the combined data set. 



SN data set Q 



M 



X 



X 2 /dof 



'gold' 0.287 + 0.026 178 1.14 

SNLS 0.265 + 0.022 114 1.00 

Combined 0.274 + 0.017 239 1.05 

We see from Table 1 and Fig. 3 that the fits to the two data sets are quite 
compatible. It is therefore reasonable to make a joint fit to the combined data 
set, which yields the results shown in the third row of Table 1: the low value 
of the x 2 /dof justifies this combination a posteriori. For the value of f2 M found 
there, the critical redshift z* above which deceleration occurs is z* ~ 0.74. 
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Fig. 4. The likelihood function x 2 an d the reduced x 2 P er degree of freedom as 
a function of ^o of the super-horizon model [6]. (i) Left figure: 'Gold' sample of 
157 supernovae [4] and SNLS sample of 71+44 supernovae [18]; (ii) right figure: 
combined 'gold' and SNLS sample of 157+71 supernovae. 

4-2.2 Super-horizon model 

For the super-horizon model of Ref. [6] we find the fit values shown in Table 2. 
Figure 4 displays the \ 2 values for the two data sets, as functions of the model 
parameter \l/£o- in this case, the consistency between the fits to the 'gold' 
and SNLS data sets is not as good as for the ACDM model. Nevertheless, 
the overall fit quality is acceptable, and the third line of Table 2 displays the 
parameter values found in a combined fit: the low value of the % 2 /dof again 
justifies this combination a posteriori. In this model, the critical redshift is 
z* = 4.6. 

Table 2 

Fits to the SHCDM model [6] parameter ^q using the 'gold' and the SNLS data 

sets and, in the third line, the result of a fit to the combined data set. 



SN data set 


*£0 


x 2 


X 2 /dof 


'gold' 


-0.79 + 0.08 


183 


1.17 


SNLS 


-0.94 + 0.09 


114 


1.00 


Combined 


-0.87 + 0.06 


245 


1.09 



The value of the critical redshift z* marking the transition from deceleration to 
acceleration may be compared in each of the above cases with that reported 
in Ref. [4]. We note that the precise value of z* depends crucially on the 
microscopic model used, in particular on the equation of state, which is in 
general z-dependent. This affects the precise functional dependence on z of 
the Hubble parameter H(z) and the deceleration parameter q(z), as seen in 
Fig. 5. The analysis of Ref. [4], which indicated a critical z* = 0.46 + 0.13, was 
based on a very simple assumed form (linear) for the z-dependent deceleration 
rate, namely q(z) = qo + q±z, which is not the case in all models, as seen in the 
right panel of Fig. 5. We therefore consider that the issue of the exact value 
of z* is a delicate point to be explored in the future. 
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Fig. 5. The Hubble parameter (left) and the deceleration parameter (right) as a 
function of the redshift z for various cosmological models: (i) red long dashed curve, 
the standard ACDM model, (ii) blue short dashed curve, the SHCDM model of 
Ref. [6], and (iii) green solid curve, Q-cosmology (off-shell) model with dark energy 
relaxing as 1/t 2 . 

4.2.3 Q- Cosmology (off-shell) Model 

We now present a detailed fit of the off-shell Q-cosmology model described 
above to the presently available astrophysical data on supernovae [4, 18], 
checking the consistency of the model and determining the crucial cosmo- 
logical parameters of this generic class of Liouville Q-cosmology models of 
Refs. [8-10]. 

We again analyze the supernovae data [4, 18], as given in terms of the distance 
modulus. We use the formula (39), with the Hubble parameter given by (36), 
to obtain di, as predicted by the model of Refs. [8, 9]. As before, we use the 
157 data points from the so-called 'gold' sample of supernovae [1, 4] observed 
by ground-based facilities and Hubble Space Telescope, as well as the more 
recent SNLS sample of 115 supernonae reported in Ref. [18]. The best-fit 
parameter values and the corresponding values for the likelihood function x 2 
are listed in Tables 3 and 4 for both datasets and their combination. 

The analysis in Table 3 has three free parameters, namely Q3, Q$ and the 
exponent 5. As in previous fits, the parameter values favoured by the 'gold' 
and SNLS supernovae samples are quite compatible, so we concentrate on 
fit to the combined data. We see that the three contributions to the energy 
density, fl2,3,s, are not very well determined individually, but the fit prefers a 
value close to 4 for the exponent S. We therefore repeated the fit, this time 
fixing (5 = 4, with the results shown in Table 4. The quality of the combined 
fit is no worse than with S left free, and the values of the different energy 
densities ^2.3,4 are compatible with those in Table 3, but with significantly 
smaller errors. The Q-cosmology curves in Figs. 1 and 2 are based on the 
best-fit values in Table 4. 9 Off-shell contributions to the 'dust' term Q-i allow 



This analysis was repeated using also the so-called 'silver' supernovae, and the 



21 



it to take negative values, so the negative signs found in the fits in Tables 3, 

4 are quite consistent with this framework. 

Table 3 

Q-cosmology model parameter values favoured by the 'gold' and the SNLS super- 
novae datasets. Q2 is determined by the other densities so that $1$ + O3 + Q2 = 1- 

SN dataset fi 3 Us 5 fi 2 X 2 X 2 / dof 

'Gold' -3.2 ±0.4 1.0 ±1.7 4.10 ±0.20 3.2 ± 1.7 177 1.14 

SNLS -2.0 ±0.8 0.4 ±0.4 4.7 ±1.0 2.6 ± 0.9 113 1.01 

Combined -3.7 ±1.1 1.3 ± 0.7 3.9 ± 0.3 3.4 ± 1.3 237 1.05 

Table 4 

Q-cosmology model parameter values favoured by the 'gold' and the SNLS super- 
novae datasets for a fixed value 5 = 4. ^2 is determined by the other densities so 
that ^4 + ^3 + ^2 = 1. 

SN dataset fi 3 fi 4 ^2 X 2 X 2 / dof 

'Gold' -3.6 ±0.8 1.2 ±0.3 3.4 ± 0.9 177 1.14 

SNLS -3.0 ±0.9 1.0 ±0.4 3.0 ± 1.0 113 1.00 

Combined -3.3 ± 0.6 1.11 ±0.25 3.2 ± 0.7 237 1.05 

We observe from Fig. 5 that, with the best-fit values of the Q-cosmology 
parameters in the model of Refs. [8, 10], the transition from a decelerating 
expansion of the Universe to an accelerating expansion would have taken place 
when z = z* ~ 0.31. This is almost compatible at the \-a level with the result 
of the analysis of Ref. [4], which indicated a critical z* ~ 0.46 ±0.13. However, 
this was based on a very simple (linear) assumed form for the deceleration 
q(z) — qo + qiz, which is not the case of our model, as also seen in Fig. 5. It 
should be possible in the future to explore observationally nonlinear forms of 
w^z), with the aim of improving the fit and reducing the uncertainty in the 
value of the critical z = z* for the transition from deceleration to acceleration. 

5 Results and Outlook 

We have shown in this paper that the two main data sets for high-redshift su- 
pernovae are compatible, at least as far as their constraints on the cosmological 
models studied in this paper are concerned. The standard ACDM model fits 
the supernova data very well, and the super-horizon dark matter model also 
fits the supernova data quite well. Both of these models are on-shell, i.e., they 
satisfy the pertinent Einstein's equations. However, off-shell cosmology models 



results were comparable to the ones presented here, confirming the robustness of 
our analysis. 
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can still be compatible with the data. As we discussed above, off-shell effects 
are important in our Q-cosmology model. Introducing an extra parameter to 
allow for these off-shell effects, we find that the Q-cosmology model may fit the 
supernovae data as well as the standard ACDM model. 

The best fits acquired by combining the 'gold' sample [4] and the SNLS [18] 
supernovae data set have the following x 2 values: 

• In the ACDM model that combines a cosmological constant with cold dark 
matter, assuming a flat Universe we find x 2 = 239 (x 2 /dof = 1.05) for 
tt M = 0.274 ± 0.017, corresponding to a value of fi A = 0.726 ± 0.017 for 
the dark energy density. This value is consistent with the one predicted in 
Ref. [4]. 

• In the super-horizon dark matter model (SHCDM) of Ref. [6], we find x 2 = 
245 (x 2 /dof = 1.09) for the best-fit case with ^ £0 = -0.87 ± 0.06 . 

• For the Q-cosmology model of Refs. [8, 10], assuming S = 4: x 2 = 237 
( x 2 /dof = 1.05): for fi 3 = -3.3 ± 0.6, fi 4 = 1-11 ± 0.25, yielding fi 2 = 
3.2 ±0.7. 

Just as this analysis was being completed, the three- year WMAP data (WMAP3) 
were released [2b]. Assuming a flat Universe, the WMAP3 data by themselves 
yield within the ACDM model the result tt M = 0.238+o;°l?. This central value 
is slightly lower than our result, but the error is larger, and the two results are 
compatible within their errors. The WMAP collaboration has combined their 
data with other CMB data, and the most precise result, obtained by combin- 
ing WMAP3 with the 2dF Galaxy Redshift Survey, is Vt M = 0.236+g;^, which 
differs from our combined result by about 2a. On the other hand, combining 
WMAP3 with large-scale structure data gives larger central values of Qu- 
The smallest errors are those obtained by combining WMAP3 with the Sloan 
Digital Sky Survey LRG data set: fi M = 0.267+H^ and with the CFHTLS 
lensing data: Hm = 0.299to;o25- These results straddle our combined fit to the 
supernova data. 

The data on baryon acoustic oscillations (BAO) provide complementary in- 
formation [19, 20]. Within the ACDM model and assuming a flat Universe, 
the BAO yield an independent estimate [20] Qm — 0.273 ± 0.024, which is 
completely compatible with the value we obtained from the supernova data. 
Combining the two results, we find Qm = 0.274 ± 0.014. We leave for a future 
occasion the exploration of the constraints imposed by the BAO data on the 
Q-cosmology models. When confronting the BAO data with the Q-cosmology 
model, it will be desirable also to formulate more precisely the predictions of 
this model, which is also an important subject for future work. There is one 
important issue which should be stressed at this point. The exotic matter scal- 
ing: a~ 6 with 5 ~ 4 of the Q-cosmology model at late eras, if valid at earlier 
times, would have dominated the epochs characterised by redshifts z > 10. As 
a result, the positions of the peaks of the Baryon Acoustic Oscillations would 
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have been modified, leading to stringent constraints on the models. However, 
in such a scenario the phenomenological observables currently used in analyses 
of the BAO [20] would be inapplicable to Q-cosmology models. We hope to 
come back to this issue in a future publication. 
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